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ABSTRACT 


The  objective  of  this  study  is  to  analyze  the  behavior  of 
a  typical  ship  structural  member  exposed  to  fire. 

Included  in  this  »/jork  are  chapters  pertaining  to  the 
modelling  of  the  fire  process  and  the  structural  member  to  be 
analyzed.  After  determining  the  expected  temperature  exposure 
of  a  shipboard  structural  member,  a  heat  transfer  analysis-  is 
carried  out  and  the  resulting  load  is  subsequently  used  in  the 
failure  anal ys i s . 


The  major  thrust  of  this  study  is  the  analysis  of  the 
failure  mechanism  of  the  collapse  of  an  I  beam.  Two  methods 
of  analysis  developed  by  Kecman,  and  Hayduk  and  Wierzbicki, 
are  summarized  and  appl  i ed  to  the  case  of  a  fire  load.  The 
reduced  moments  due  to  fire  exposure  are  determined  and 
finally  recommendations  are  offered  for  the  use  of  this 
information  in  ship  structural  design. 


Thesis  Supervisor:   Tomasz  Wierzbicki 
Title:   Professor  of 

Applied  Mechanics 
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Chapter  1 


iNTRnnnnTTnN 


Traditionally,  Marine  Engineers  have  handled  the  effects 
of  fire  on  shipboard  structures  through  the  use  of  damage 
control  procedures.  Ideally,  a  structure  should  be  designed 
to  withstand  anticipated  fire  loads.  The  objective  of  this 
work  is  to  apply  the  methods  of  post  critical  analysis  of 
thin-walled  structures  developed  by  Mechanical  Engineers  to 
analyzing  fire  loads  on  board  ocean  engineering  structures, 
specifically  ship  structures. 

A  ship  structure  is  a  highly  redundant  structure.  If  the 
load  bearing  capacity  of  one  member  is  diminished,  then  the 
surrounding  structural  members  carry  more  of  the  applied 
load.  Shipboard  fires  are  characteristically  highly  localized 
phenomenon  which  may  effect  only  a  small  percentage  of  the 
entire  structure.  By  examining  the  effect  of  fire  on  a  single 
structural  member,  one  will  gain  insight  into  the  behavior  of 
the  entire  structure.  Ultimately,  ship  structures  may  be 
designed  to  withstand  the  effects  of  typical  shipboard  fires. 

It  is  the  intent  of  this  study  to  examine  the  effect  of  a 
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■fire  on  a  typical  shipboard  structure,  namely  the  I  beam.   The 

major  thrust  will  be  the  analysis  of  the   -failure  mechanism  of 

the  beam.   Two  methods  will  be  introduced  which  analyze   the 

post   yield   behavior   o-f   an   I   beam.    Figure  1.1    shows 

graphically    how   a   -fire   load   diminishes   the  mechanical 
properties  o-f  a  beam. 

The  material  yields   upon   reaching  the  maximum  moment  in 
Figure  1.1.  Consider  the  -following  relation: 

Sti-f-fness  =  Moment/Curvature. 

Beyond   the   yield  point,  sti-f-fness  is  less  than  or   equal   to 

zero.   Normally  a  structural  member   remains   in   the   elastic 

range,  and  the  moment  remains  less  than  M  ,   the  -fully  plastic 

o 

yield   moment.    However,   through   exposure   to   fire,    the 

sti-f-fness  will  decrease  so  rapidly  that  the  moment   will   -fall 

below  a  sa-fe  range. 

The  ultimate  objective  o-f  this  study  is  to  show  how  much 
the  strength  o-f  an  I  beam  is  diminished  when  exposed  to  -fire. 
To  this  end,  the  structure  and  the  fire  process  will  be 
modelled;  and  the  effect  of  temperature  on  the  mechanical 
properties  of  a  selected  material  will  be  given.  After  these 
models  have  been  developed,  a  heat  transfer  analysis  will  be 
performed.  As  a  result  of  this  analysis,  it  will  be  shown 
that  bending  collapse   is   the  primary  mechanism  of  structural 
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■failure.  Next,  the  failure  mechanism  analysis  methods  will  be 
■fully  developed.  Finally,  recommendations  will  be  o-f-fered  -for 
applying  the  results  o-f  this  work  to  the  design  o-f  an  entire 
sh  i  p  struc  ture  . 
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Chapter  2 
MDDF1  DF  STRIinTIIRF 


This  chapter  will  develop  the  motivation  towards  choosing 
an  I  beam  as  a  typical  shipboard  structural  element. 
Additionally,  a  material  will  be  selected;  and  the  temperature 
dependent  properties  o-f  that  material  will  be  given. 


Z-.l_BajiJig£:iDUBd 

The  principal  aspects  o-f  ship  structural  design  most 
commonly  employed  by  Marine  Engineers  involves  the  application 
o-f  rules  established  by  Societies  such  as  the  American  Bureau 
o-f  Shipping.  These  rules  are  based  upon  the  primary  loads 
experienced  by  a  ship  during  its  lifetime.  These  loads 
include:  still  water  bending  moments  and  shearing  -forces;  wave 
induced  moments  and  shearing  -forces;  impulse  loads;  thermal 
loads;  and  dynamic  loads. 

These  loads  which  are  unique  to  ocean  going   vessels   are 
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often  difficult  to  accurately  predict.  There-fore,  ship 
structural  design  is  largely  a  process  involving  guesses 
derived  -from  intuition  or  empirical  data.  The  primary  stress, 
strain  and  displacement  behavior  of  the  ship  structure  are  due 
to  displacements  limited  to  those  parallel  to  the  unloaded 
midsurface  of  the  plate.  Thus,  traditionally  hull  strength  is 
determined  through  an  analysis  of  the  midships  section. 
Figure  2.1  illustrates  a  typical  midships  section. 


2*  2— Load— In  AD3iru-s.sl.Qn 

Today,  Ship  Structural  Engineers  have  adopted  a  model 
which  is  motivated  by  considering  the  loads  transmitted  to  the 
structure.  Loads  transmitted  to  a  stiffened  plate  (by 
slamming,  etc.)  are  in  turn  transmitted  to  the  stiffeners  and 
from  them  to  the  point  supports,  and  so  on.  Eventually  a  load 
transmitting  grillage  is  derived.  Next,  boundary  conditions 
will  be  used  to  define  if  the  spans  are  treated  as  a  series  of 
isolated  points  or  continuous  beams.  In  most  cases  the 
boundary  points  will  be  those  which  can  be  described  as  simple 
or  clamped  supports. 

Figure  2.2  shows  examples  of  typical  members  o-f  a  ship's 
structure.   For  the  purposes  of  this  study,   an  I  beam  will  be 
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analyzed  as  its   behavior  is  considered  representative  of  the 

open  beams  in  Figure  2.2.  Figure   2.3  gives  the  dimensions  o-f 

the  beam  which  herea-fter  will   be   used   as   the   model  to  be 
anal yzed. 

The  -first   moment   of  area  o-f  the  I  beam  will  be  required 
i  n  the  anal ysi  s . 

b2t 

/  zdA  =  t  (b  -  t_)  (a  -  t  )  +  2£  (2.1) 

H  r      r       w     4 

A  =  Sec  t  i  on  Area 

t   =  Thickness  o-f  the  -flanges 

-f 

t   =  Thickness  o-f  the  web 
w 

a  =  Width  of  flanges 

b  =  He  i  gh  t  of  web 

2 
Assume  t   is  negligible.   Hence,  the  first  moment   of  area  may 


be  approximated  by: 


2 

r     ^  w  <2.2) 

_/  zdA  -   abt^  +   - 
A  f     4 


2^2_3&le.c±xon_o£_td-a±^clal 

Several  materials  have  been  used  throughout  shipbuilding 
history  starting  with  wood,  and  progressing  to  iron  and 
steel.   Today  several  grades  of  steel  and   even   titanium   are 
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used  in  the  shipbuilding  industry.  One  of  the  most  common 
grades  o-f  steels  used  is  ABS  Grade  A.  Material  Scientists, 
however,  have  not  tested  ABS  Grade  A  in  order  to  determine  the 
temperature  dependent  material  properties.  Ordinary 
structural  steels  such  as  ASTM-A36,  have  been  widely  used  as  a 
test  specimen.  Table  2.1  compares  the  properties  o-f  the  two 
stee 1 s. 


Table  2.1 

St 

eel  Properti 

ASTM-A36 

0.23 

0.04 

0.05 

36.00 

i  es 

ABS  Grade-A 

Carbon,  max  V,  0.23  0.26 

Phosphorous,  max  '/.  0.04                       0.04 

Sulpher,  max  V.  0.05  0.0  4 

Yield  Stress,  ksi  36.00  34.00  (min> 


It  is  considered  a  conservative  assumption  that  an 
analysis  -for  ASTM-A36  can  be  used  without  sacr  i  f  i  c  i  ng 
accuracy ,  <  ref  .  1 ) . 

One  of  the  major  contributions  to  the  loss  of  strength  is 
the  temperature  dependence  of  the  mechanical  properties  of  the 
materials  under  consideration.  Several  researchers  have 
predicted  the  deterioration  of  steel  structural  elements  at 
elevated  temperatures.  A  search  of  the  literature  reveals 
good  agreement  among  those  whom  have  studied  the  temperature 
dependence  of  material  properties. 
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Harmathy  and  Stanzak  <ref.  2)  have  carried  out  several 
experiments  in  order  to  provide  to  Design  Engineers,  primarily 
to  those  concerned  with  -fire  endurance  of  structural  elements, 
the  elevated  temperature  mechanical  properties  of  common 
commercial  steel  products.  The  temperature  dependent 
expression  -for  Young's  Modulus  o-f  Elasticity,  E,  <psi)  for 
ASTM-A36  steel  is  given  here: 

E  =  30*106  -  9.3T"  (2.3) 

o 
T  =  Temperature  of  heat  exposure,   Fahrenheit 

Steel  loses  strength  as  it  is  heated.   Figure  2.4   illustrates 

the  temperature  dependent  tensile   strength  of  ASTM-A36  steel. 

o         o 
Notice   that   between  1000   and  1200  F   the   tensile   strength 

falls  below  the  usual  design  stress. 

For  the  purposes  of  this  study,  linear  dependence  of 
temperature  is  assumed.  The  following  equation  will  be  used 
in  this  anal.ys  i  s  . 

O   =  -20. 5T  +  40,000.0  (2.4) 
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A  survey  of  the  literature  reveals  several  methods  of 
characterizing  a  -fire.  Harmathy  offers  a  numerical  technique 
which  is  summarized  in  Section  3.1.  Most  Fire  Engineers  use  a 
time  temperature  curve  to  define  the  characteristics  of  a 
fire.  The  development  of  an  acceptable  time  temperature  curve 
is  given  in  Section  3.2. 


3^1_£.ha.cax:l.a.El.sl_Lc.s._.Qj£._a—  Shlfib&and-ElnB. 

The  heat  exposure  on  a  shipboard  structure  resulting  from 
a  fire  is  a  function  of  several  identifiable  and  incidental 
variables,  (ref.  3).  The  three  primary  factors  contributing  to 
the  severity  of  fully  developed  fires  is  the  flow  rate  of  air 
into  a  burning  compartment,  (compartment  ventilation);  the 
total  fire  load;  and  the  thermal  characteristics  of  the 
compartment  boundaries. 
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Numerous  theoretical  studies  have  determined 
time-temperature  curves  of  typical  fires  based  on  these 
factors,  (ref  3).  Even  with  simplifying  assumptions,  these 
factors  are  only  calculable  with  a  computer.  Harmathy 
introduces  three  easily  calculable  -fire  severity  parameters 
which  can  be  assumed  to  be  accurate  predictions  of  the 
contributing  factors  to  a  -fully  developed  fire.  These 
parameters  are  defined  as  follows: 

1.  The  duration  of  the  fully  developed  fire,  t; 

2.  The  overall  penetration  flux,  q,  that  is  the  heat  flux 
absorbed  by  the  compartment  boundaries,  averaged 
spatially  over  the  boundary  surfaces,  and; 

3.  The  average  temperature  of  the  compartment  gases, 

(average  "fire  temperature"),  T  ,  averaged  spatially  over 

g 

the  period  of  full  development. 

It  is  useful  at  this  time  to  introduce  the  three  periods 
of  freely  burning  compartment  fires.  See  Figure  3.1.  The 
focus  of  this  study  is  on  the  period  of  full  development, 
since  80-90  percent  of  the  fuel  energy  is  released  during  this 
per  i  od . 

A  compartment  fire  is  a  fire  which  is  enclosed  in  a  shell 
of  arbitrary  shape.  Typically  on  board  a  ship,  ventilation  is 
secured   to   the   compartment   with    the    fire,    and    free 
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communication  to  the  surrounding  spaces  is  not  allowed.  Fire 
load  is  a  numerical  expression  which  is  a  -function  of  the 
available  -fuel  in  the  compartment.  The  units  ot"  specific  fire 
load  is  pounds  of  wood  per  square  foot,  <ref.  1). 

To  compute  the  fire  load  of  a  particular  compartment,  the 
products  of  the  weights  in  pounds  of  the  material,  and  the 
heat  capacity  are  summed  for  all  the  combustible  materials  in 
the  compartment.  The  resulting  expression,  which  has  units  of 
BTU's  must  be  divided  by  8000  BTU/lb.  This  result  will  provide 
the  same  available  energy  as  the  equivalent  weight  of  wood. 

Equation  <3.1)  gives  the  equivalent  weight  of  wood  in 
pounds  divided  by  the  deck  area  of  the  compartment. 


F  =  E  W . H . /H  A  <  3  .  1  > 

11   w  F 


2 
F  =  Specific  fire  load,  lb/ft 

Ul   =  Weight  of  individual  materials,  lb 
i 

H   =  Heat  capacity  of  individual  materials,  BTU/lb 
i 

H   =  Heat  capacity  of  wood,  BTU/lb 
w 

2 
A   =  Area  of  deck ,  ft 

F 


Table   3.1   gives  typical  fire  load   values   for   various 
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shipboard  compartments,  (re-f.  4). 


Table  3.1 

Spec  i -fie  Fire  Loading  o-f  Shipboard  Spaces 

lbs. /square  -feet 
Control  Spaces 

Ulhee  1  house/Char  troom  1.5 

Fire  Control  Stations  1.5 

Escape  Routes 

Com  i  dors  1  .5 

Stairway  Enclosures  1.0 

Ac c ommoda t  i  on  Spaces 
Staterooms 

Fire  Resistant  Furnishings  3.0 

Combustible  Furnishings  5.0 

Pub!  i  c  Spaces 

Fire  Resistant  Furnishings: 

Lounges,  Restaurants,  etc.  3.0 

Ferry  'Jesse  Is  1.5 

Combustible  Furnishings  5.0 

Sanitary  Spaces  0.4 

Serv  i  ce  Spaces 

Gal  1 eys  10.0 
Pantries  (no  heating  appliances)  4.0 
Food  Concession  (-ferry  vessels- 
no  combustible  stowage)  1.5 
Workshops  10.0 
Storerooms  -  Combustibles  10.0 
Storerooms  -  Cleaning  gear  only  3.0 
Ship's  Laundry  10.0 
Laundry  -  Private  Use  1.5 

Main  Machinery  and  Cargo  Spaces  10.0 

Auxiliary  Machinery  Spaces  5.0 

Tanks  and  Voids  0.0 
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It  has  been  determined  by  several  researchers  that  the 
■flow  rate  of  air  into  a  burning  compartment,  (ventilation)  is 
just  as  critical  to  the  severity  of  the  fire  as  the  fire 
load.  The  determination  of  compartment  ventilation  depends  on 
how  the  rate  of  entry  of  air  compares  with  a  critical  value. 
The  rate  of  entry  of  air  into  a  compartment  can  be  expressed 
as: 


U   =  C,A  h 
a    1  w  w 


1/2 


<3.2> 


U   =  Flow  rate  of  air,  1  b/hr 


A   =  Area  of  window,  ft 
w 

h   =  Height  of  window,  ft 


w 


5/; 


C   =  240  lb/ft    -hr 

1 


The  critical  value  is  defined  as: 


(U_)__  -  c24>faf 


a  cr 


(3.3) 


C   =23.1  lb/ft  -hr 
2 

o   =  Material  constant  which  depends  on  the  furnishings  of 
the  compartment.   For  typical  living  spaces, 
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0  .55  <  9    <  0 .90. 

The  air  flow  factor  is  expressed  as  -follows: 

u 


a  cr 


Thus,  a  -fire  is  ventilation  controlled  if  £  <  1.0;  and  a 
fire  is  fuel  surface  controlled  if  £  >=  1.0. 

The  average  rate  of  burning  of  combustibles,  R,  has  been 
determined  during  the  period  of  full  development  based  upon 
empirical  data.   For  £  >=  1.0 


R  =  C3$FAF  (3.5) 


2 
C  =4.57  lb/ft  -hr 
3 


For  l    <  1.0 


R  =  0.163U 

3. 


(3.6) 


The  mass   flow   rate   of   gaseous   products   out   of   the 
compar  tmen  t  is: 

U  =  R  +  U  (3.7) 

g        a 
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It  -foil  cws  that  : 


Q  =  R(0.9326H   +  0.0686H  )  (3.8) 


where  : 


Q  =  average  rate  of  heat  evolution  -from  combustible 
mater  i  al s 

H   =~  5760  BTU/lb,  heat  of  combustion  of  char 

1 

H   =~  11480  BTU/lb,  heat  of  combustion  of  volatile 
2 

decomposition  products 

5  =  the  fraction  of  heat  released  by  the  combustion  of 
the  volatile  decomposition  products  inside  the 
compar  tmen  t 

6  is  calculated  based  upon  an  estimation  of   the   average 
f 1 ame  length,  1 : 


1/3 
C   =  1 .1 1  ft 
4 


I  =   ccu  1/3  (3.10) 

5  a 


1/3      1/3 
C   =  0 .366    f t(h)    /(lb) 
5 
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Then  , 


if  £  =  Hc  6  =  1  <  3 . 1  1  ) 


if   I  >  H  6  =    (H  /£)1/3  <3.12) 

c  c 


where    H       is    the    height    of    the    compartment. 

c 


The    duration    of    a    fully    developed    fire,     t,     is      calculated 
as    f ol 1 ows : 


for  £  =  1  j  =  C  /<$>  (3.13) 

6 


2 
C      =    0.205    ft    <hr)/lb 
6 

f or  £   <    l  T  =  c  /^  (3.14) 

6 


The  other  two  parameters  may  be  obtained  numerically  from 
the  following  equation: 


1  A  A 

q  =  -  [p  -  U  c  (CT   -  T  )  -  0-^(T    -  T  4)] 

At     ggg   o     3  g    o 


(3.15) 


where  : 


A   =  Total  area  for  the  compartment 

t 

Q  =  Average  rate  of  heat  evolution  (BTU/hr) 

o 
C   =  Specific  heat  of  compartment  gases  <BTU/lb-  R) 
9 


-  2?  - 


o 
T   =  Temperature  of  combustion  gases,  R 
9 

o 
T   =  Temperature  at  time  =  0,  R 
o 

-8        2      o   4 
0  -  Stefan-Bol tzman  constant,  0.171*10    BTU/ft  -hr-C  R] 


£  =  Empirical  factor,  1.05,  d i mens i on  1  ess 


Tg  -  ts_  +  (Tq  +  a(2TKA)vW/4 


where  : 


>7  =  Empirical  -factor,  0.9,  d  i  mens  i  on  1  ess 

2 
k    =  Average  thermal  diffusivity,  ft  /hr 

o 
k  =  Average  thermal  conductivity,  BTU/ft-hr-  R 


Figures  3.2  and  3.3  show  the  dependence   of   the   overall 

penetration   flux,    and    average    temperature    of    gases, 

respectively,   on   the  air  flow  factor.    These   figures   were 

2 
determined  empirically   for  a  specific  surface  of  0.63  ft  /lb. 

Curves  are  given  for  three   values   of  the  specific  fire  load: 

2 

3.1 ,  6.1 ,  and  12.3  1 b/f t  . 

Figure   3.2   shows  clearly  that  the   maximum   penetration 
flux,  and   fully   developed  fire  duration  occurs  at  lt   =  1.  The 
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specific  -fire  loading  appears  to  have  little  effect  on  the 
maximum  -flux,  and  the  duration  is  usually  quite  short.  It  may 
also  be  concluded  from  Figure  3.3  that  the  hottest  -fires  occur 
in  ventilation  controlled  compartments.  It  is  possible  at 
this  time  to  draw  some  conclusions  about  shipboard  fires. 
Below  deck's  fires  which  do  not  communicate  freely  with  the 
atmosphere  are  ventilation  controlled.  These  fires  are  more 
destructive  as  the  duration  is  greater,  the  average 
temperature  is  greater,  and  the  effective  heat  flux  is  greater 
than  for  an  above  decks  fire. 

McCarthy    demonstrated   the   usefulness   of    Harmathy's 

technique  for  predicting   fire   severity  of  a  typical  berthing 

compartment.   Figures   3.2   and   3.3  can  be  used  for  a  typical 

berthing   compartment,   and  in  accordance  with   Table   3.1,   a 

2 
specific  fire   load   of   3.0  lb/ft   was  used.   He  assumed  that 

the  ventilation  is  secured,  and   that   one  water  tight  door  is 

open  for  access  to  -firefighters. 

2 
A   =  15.0  ft 
w 

h   =  10  ft 
w 

U   =  11334.2  lb/hr 
a 

2 
A   =  150  ft 

F 

<P    =  0  .75 
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CU  ]    =  9483.75 
a  cr 

£  -  1 .20 


Hence,  the   fire  is  fuel  surface  controlled.   Notice  that 

as   A    increase  the  fire  approaches   ventilation   controlled. 

F 

From  Figures  3.2,  and  3.3: 

3 
t  =  1.15*10   sec  =  19.2  min 

2 
q  =  6344  BTU/hr-ft 

o 

T   =  1256  F 
9 


It  should  be  appreciated  from  this  discussion  that  fire 
is  indeed  a  complex  process  with  several  unquan t i f i abl e 
variables.  The  term  destructive  potential  was  introduced  in 
reference  1.  The  destructive  potential  of  a  fire  may  be 
considered  a  function  of  the  maximum  temperature  obtained,  the 
duration  of  the  fire,  the  duration  of  time  at  maximum 
temperature,  and  the  quantity  of  energy  transmitted  from  the 
fire.  Most  Fire  Engineers  use  a  time  temperature  curve  as  a 
definition  of  heat  exposure.  For  the  purpose  of  this  study, 
the  use  of  an  accepted  time-temperature  curve  will  be 
considered  an  adequate  expression  of  the  fire  process. 
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Traditionally,  -fire  tests  have  been  performed  based  upon 
prior  established  standard  time  temperature  curves.  The  most 
widely  used  fire  test  in  this  country  is  the  ASTM-E11?  which 
was  developed  by  the  American  National  Standards  Institute, 
(re-f.  5).  This  method  was  established  as  a  result  of  the 
necessity  to  standardize  the  structural  performance  of 
building  construction. 

The  basic  procedure  of  any  fire   test   is   to   place   the 

element  to  be  tested  into  a  gas   furnace  in  which  the  standard 

time  temperature   curve  is  duplicated.   The  gas  temperature  is 

closely   monitored    through    the    use    of    thermocouples. 

Additionally,  deflection,  axial  thrust,  and  collapse   behavior 

are  measured.   The  Ell?  time  temperature  curve  initially  rises 

o 
rapidly  reaching  1000  F  at  5  minutes  into   the   test.    Beyond 

this  point,  the  temperature  is  increased  more  slowly   so   that 

o 
at  the  end  of  one  hour  the   furnace  temperature  is  1700  F.  The 

o 
maximum  test  temperature,  2300  F  is  reached   at   5   hours   and 

held  constant  thereafter. 

Another  test   method   commonly   used   by   Ship  Structural 
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Engineers  is  the  International  Convention  -for  the  Safety  of 
Life  at  Sea  (SOLAS)  A60  test.  This  test  is  used  for  testing 
non-structural  bulkheads  from  berthing  and  office 
compartments.  The  specimen  observed  is  not  loaded  but 
strained  at  the  edges.  Other  non-structural  bulkhead  tests 
are  SOLAS  A30  ,  A15,  AO  ,  B15,  BO,  but  the  requirements  for 
these  tests  are  not  as  severe  as  for  the  A60  test,  <ref.  1). 

For  the  purposes  of  this  study,  reference  6  supports   the 

o 
use  of  a  step  change   to   1700  F  vice  the  gradually  increasing 

temperatures  of  ASTM-E11?  or  the  SOLAS  A60  curves.   Figure  3.4 

gives  the  ASTM-E119   curve,   the   SOLAS   A60   curve,   and   the 

assumed   time   temperature    relation    to   be   used   in   this 

anal ys  i  s . 
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Chapter  4 
HFAT  TRANSFFR  ANA1YSTS 


The   intention   of   this   analysis   is   to  determine   the 

relation  between  temperature  and  the  distance  along  the  length 

of  the  beam.   A  beam  exposed  to  a  fire  load  will   be   modelled 

as  a  fin  with  a  temperature  gradient.   The   temperature  at  the 

exposed  face  is  T  ;  the  temperature  at  the  end  of  the   fin   is 

b 

T  ,  the  ambient  temperature.    The   geometry   of   the   fin   is 

a 

defined  in  Figure  4.1 


4-.J Ccinduc,±lon_d£a±_Ic:ans££ic. 


The  following  assumptions  are  used: 


-  T   =  Temperature  at  x  =  0 

b 

-  T   =  Temperature  at  x  =  L 

a 

-  The  heat  transfer  is  steady  state. 
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-  k  =  Uniform  conductivity  o+"  the  material, 


o 
BTU/E  in-sec-  F3 


-  h  =  Uniform  surface  coefficient  of  heat  transfer, 


BTU/Chr-f t  -  F3  . 

Additionally,  this  study  will  use  the  assumption  that  the 
rod  is  thin,  which  implies  that  the  temperature  in  the  thin 
direction,  perpendicular  to  x,  is  uni-form.  Therefore,  the 
surface  temperature  equals  the  centerline  temperature,  (re-f. 
7).  Applying  an  energy  balance  on  an  elemental  section  *x , 
yields  a  governing  differential  equation.  In  steady  state, 
with  no  internal  heat  generation: 

Expanding  the  first  term  of  the   right   hand   side   gives   the 
following  equation: 

Sc   ~    Sc+Ax   =   %   "    (qx   +    8x2lAx)  (4.2) 

-8qx 
^        qx+Ax        IbT^*  <4.3) 

%   "   qx+Ax   =   "—  C-kAyAz&Ax  <4.4) 

dx  ox 

Sc  -   Sc+Ax  =i-(k&^)Ax  <4.5) 

dx  dx 

Uhere    S   =    Surface    area   =   A/Az.    Conduction    heat    transfer    across 
.A  A    is    given     in    the    following 
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equat  i  on  : 


h  =  ^  (4.6) 


T-T 
a 


q.       =    h(T    -    T    )A  (4.7) 

AA  a 


q^   =    hPAx(T   -    T    )  (4.8) 


Where 


P   =    Perimeter    =    20y    +    _\z) 
T   =    Temperature    of    fin    at    x 

Substituting    (4.5)    and    (4.8)     into    (4.1)    yields: 


d_(ks^I)Ax  -   hPAx(T  -  T   )    =  0  (4.?) 

dx        dx  a 


Finally,  the  governing  differential  equation  for  6  =    CT  - 


T  3  is: 
a 


^    h£  ft    n 

,2    kS  (4.10) 

dx^ 


Where  S,  k,  and  T   are  constant. 

a 


The  differential   equation   for   one   dimensional   steady 
state  heat  flow  with  one  internal  heat  generation  source  is: 

Sl±  ~   ~  _±  (4.11) 

dx2     k 

Note  that   the   fin   problem  is  similar  to  the  single  internal 
heat 
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generation  source  problem  with  a  negative  heat  source: 


W  =  -  —  (4.12) 

i      S 


The  method  of  solution  will  follow  a  method  provided  in 
reference  8.  Assume  the  solution  is  of  the  form  given  in 
equation  (4.13): 

a  ~      Bx       -Bx 

9  =  C  e   +  re 

1       2  (4.13) 

Where  : 

B  =  [hP/kS]1/2  (4.14) 

The  boundary  conditions  for  a  fin   of   finite   length,   L,   as 
defined  bx  Figure  4.1  are: 

@  x  =  0;    6  =  t.  -  T  =  0.  (4.15) 

b    a    b 


@  x  =  L;     6  =  T   -T   =0 
a     a 


(4.16) 


Substituting  (4.15)  and  (4.16)  into  (4.13)  yields 


eb  =  c1  +  c2  (4.i7) 


BL       -BL  /a     \  o n 

0  =  C  e    +  c2e  <*.!»> 
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Sol  v  i  rig    -for    C    : 
1 

-BL 


-o   e 
C     =       b  (4.1?) 

BL      -BL 
e      -e 


C      -follows    from    (4.17).    Finally    the    solution     is: 
2 

9  eB(L'x)    +  e-B(L~x)  (4.20) 

6,  BL  -BL 

b  e        -    e 

At  this  point,  it  is  necessary   to  solve  -for  B.  For  the  I 
beam  in  Figure  2.3,  the  perimeter  to  surface  area  ratio  is: 

P    2(b  +  2a  -  t  ) 

S  =  w  (4.21) 

bt   -  2t_t   +  2t,a 

w      f  w      f 

Figure  4.2  gives  the  relation  between  thermal  conductivity  and 

temperature  for  plain  carbon  steel,  (ref.  ?).  Within  the  range 

o  o 

of  interest,  70  F  <=  T  <=   2000  F,   the   curve   is   considered 

approximately  constant.   For   the   purposes   of  this  analysis, 

o 
the  thermal  conductivity  at  1000  F  wi 1 1  be  assumed. 

k  =  0.0005  BTU/in  sec  °F  (4.22) 


The  condition  of  natural  convection  in  a  gas   is   assumed 
for  this  analysis,  (ref.  7): 

h  =  2.0   BTU/hr  ft2  °F  (4.23) 

Therefore , 

B  =  (P/10.8S)172  (4.24) 

where,  P/S  is  given  in  equation  (4.21) 
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The  validity  of  thin  fin   theory   must   now   be   checked. 
According  to  reference  7,  the  thin  fin  theory  is  valid  if: 

k  > 

hb=3  <4-25> 

Substituting  equations   <4.22)   and   (4.23)   into   C4.29),  the 
maximum  value  for  b  is: 

b  =  3.6  ft 

<4.26) 

This  is  a  reasonable  criteria  for  this  study.   Therefore   thin 
t*  i  n  theory  is  valid. 


4^3_S&sul±s 

It  i s  now  possible  to  draw  some  conclusions  about  the 
heat  transfer  through  the  beam.  A  Microsoft  BASIC  computer 
program  was  written  in  order  to  examine  several  geometric 
configurations  of  the  beam.  The  source  code  for  the  program 
is  in  Appendix  A.  The  results  are  given  graphically  in  Figure 
4.3.  The  dimensions  of  the  beams  examined  are  in  Figure  4.4. 

.The  following  conclusion  may  be  drawn  from  Figure  4.3.  As 
P/S,  and  L  increase,  the  effect  of  the  fire  becomes  more 
localized.   With  a  typical  shipboard   I   beam,   the   ratio   of 
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perimeter  to  sur-face  area  is  such  that  the  -fire  load  may  be 
treated  as  a  point  load.  That  is  to  say  that  the  heat 
trans-fer  through  an  I  beam  remains  highly  localized  in  the 
vicinity  o-f  the  -fire.  All  e-f-fects  o-f  the  -fire  virtually 
disappear  at  a  distance  which  is  a  quarter  o-f  the  length  o-f 
the  beam  away  -from  the  -fire. 
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Chapter  5 
ANAIYSIS  f)F  THF  FAI1  I1RF  MFHHANTSM 


As  was  stated  in  the  preceding  chapter,  an  I  beam  exposed 
to  a  typical  shipboard  -fire  will  experience  the  same  -failure 
mechanisms  as  a  beam  under  a  point  load.  Recently  there  has 
been  increasing  interest  in  this  type  of  -failure  due  to 
mechanical  loads.  It  is  the  intention  o-f  this  study  to 
approximate  the  shipboard  -fire  load  as  an  additional 
mechanical  load,  and  to  use  two  o-f  the  recently  developed 
techniques  -for  analyzing  the  -failure  mechanism.  Additionally, 
this  chapter  will  include  summaries  o-f  literature  which 
•focuses  upon  bending  collapse  of  thin  walled  structures. 

Kecman  Cref.  10)  studied  the  bending  collapse  behavior  of 
rectangular  and  square  section  tubes  both  theoretically  and 
experimentally.  He  used  a  limit  analysis  technique  and  a  set 
of  formulae  to  derive  the  relation  between  the  hinge  moment 
and  the  associated  angle  of  rotation.  Very  good  agreement  was 
found  between  his  theoretical  predictions  and  the  experimental 
results.  His  theoretical  predictions  were  in  the  form  of  a 
numerical  solution.  However,  Hayduk  and  Wierzbicki  (ref.  11), 
gave   a   closed  form  analytical  solution.    In   reference   11, 
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lower  and  upper  bound  solutions  -for  the  mean  crushing  strength 
of  cruci-Forms  were  obtained  by  considering  bending  and 
extensional  modes  of  de-formation. 


5-^1— Ka-cm.^  n-l^_tle±.hod 

Kecman  developed  a  numerical  scheme  to  analyze  the 
bending  collapse  of  rectangular  and  square  section  tubes, 
<ref.  10).  The  primary  failure  mode  is  local  buckling  and 
hinge  development.  He  considered  only  uniaxial  bending 
col  1  apse . 


5.1.1  Maximum  Bending  Strength 

Kecman  referred  to  a  study  which  examined  the  bending 
strength  of  a  box  shaped  section  with  the  webs  in  tension  and 
the  top  flange  in  compression.  The  effect  of  the  webs  was 
found  to  be  negligible.  Therefore,  the  critical  stress  of  an 
I  beam  may  be  assumed  to  be  the  same  as  the  critical  stress  of 
a  box  shape  beam. 


The  critical  stress,  a       ,  of  the  section  was  derived  in  a 

cr 
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•form : 


2 
7T  E 


cr 


12(l-v  ) 


r(5.23  +  0.1 6£)  (-) 
b  a 


<5.1) 


a  =  f 1 ange  wi  dth 

b  =  web  he  i  gh  t 

t  =  web,  and  flange  thickness 

If  the  critical  stress  is  less  than  the   yield   stress,  0      of 

o 

the  material,  then  the   -flange   buckles.    The  effective  width 

becomes: 


a   =  a(0.7-H^  +  0.3) 
e        c 
e 


<5.2) 


a       =  stress  at  the  edge  of  the  flange 

e 

If  the   critical   stress  is  greater  than  or  equal  to  the  yield 

stress,  then   a    =  a.   For  the  purposes  of  this  study,  assume 


that  collapsing  occurs   when  a         =  o    .   It   follows  that  the 

e       o 

maximum   moment,    M    ,    may    be    calculated    using    <5.1)       and      (5.2). 

m 


If    a  <=    a    : 

cr  o 


„    2a   +   b   +   a    (3-  +   2) 

M     =  a  tb2* Z* 

m  o  3  (a   +  b) 


(5. 3a) 


I  f    o         >=    3  <r  : 

cr  o 


Mm  =  Mo  =  ^t^^  ~  t)    +  0.5(b  -   2t)    ] 


<5.3b) 
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where  M   =  the  -fully  plastic  yield  moment  of  the  section.   If 
o 

o       <    o         <  3  a   ,  a  moment,  M  ',  is  defined: 
o     cr      o  o 


M  '  =  a  tb(a  +  — ) 
o     o       3 


With  the  resulting  maximum  moment: 


<5.3c) 


\  =  V  +  <M  -  M  •)  Ccr  "  Q° 


(5.3d) 


5.1.2  Actual  Bending  Collapse  Mechanism 

A  detailed  analysis  of  the  hinge  deformation  pattern 
offers  insight  into  the  failure  mechanism.  The  deformed 
pattern  is  given  in  Figure  5.1. 

The  following  characteristics  of  the  deformation 
mechanism  were  given  in  reference  10: 

(a)  Most  of  the  plastic  deformation  is  concentrated  along 
the  stationary  yield  lines:  JF ,  HI,  EB,  GB,  FC,  HC ,  CR, 
BA,  GK,  and  EL. 

<b)  Travelling  yield  lines  of  rolling  deformation  were 
also  observed  along  GA ,  AE ,  KA ,  and  LA  while  the  distance 
BA  changes  during  hinge  rotation. 

(c)  The  corner  angle  along  KL  remains  almost  constant. 
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<d)  The  deformation  of  the  tension  flange  includes  the 
travel  ling  yield  1  i  nes  NP,  and  MQ,  approx  imate 1 y 
corresponding  to  the  intersection  of  the  undeformed 
flange  planes  and  the  current  cylindrical  surface  TU.  See 
Figures  5.1,  and  5.2. 

<e)  In-plane  deformation  of  the  walls  was  also  observed, 
particularly  along  the  rolling  lines  6A,  AE,  BA,  HC,  CF , 
IR,  and  JR.  This  deformation  is  essential  to  maintain 
h  i  nge  k  i  nemat  i  cs. 

<f)  The  rolling  radius,  r,  varies  slightly  along  6A, 
measured  by  cutting  at  right  angles  to  GA  and  its 
estimated  value  was  approximately  proportional  to  the 
h  i  nge  1 ength . 

r  =~  < 0.03-0. 05) *KL 

The  radius  reduces  during  collapse,  so  that  rolling 
becomes  more  and  more  difficult,  until  the  mechanism  jams 
and  starts  a  new  phase  of  collapse. 

<g)  The  rolling  radius  along  KA,  and  LA  varies  -from 
infinity  at  K  and  L  to  a  value  close  to  r  in  <f)  above. 

Four   distinct   phases   in   the   hinge   development   were 
i  den  t  i  f  i  ed : 


<1>  Initial  phase  with  a  bulge  in  the  web,  A  ,  in  Figure 

1 
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<2)    Second    phase    with    rolling   as    described   by    the 

cross-section    shape    A       in    Figure    5.3. 

2 

<3)  Third  phase  indicated  by  a  jamming  of  the  rolling  and 

described  by  the  cross  section  B  A  D  in  Figures. 3. 

3  3 

<4)  Final  phase  initiated  by  the  contact  of  the  two 

buckled  halves  o-f  the  compression  -flange.   There  is  no 

clear  transition  from  one  phase  to  another. 


5.1.3  Theoretical  Bending  Collapse  Mechanism 

Since  structures  most  commonly  experience  hinge   collapse 

described  by  the   second   phase  which  includes  rotation  angles 

in  the  range,  5-10  <=  a   <=  25-35.  Kecman  based  his  theoretical 
analysis  on  the  following  assumptions: 

-  The  web  and  the  flanges  of  the  section  deform  along  the 
concentrated  yield  lines  only; 

-  The  web  and  flanges  are  incompressible  and 
i  nex  tens  i  bl  e  ; 

-  The  structural  continuity  is  maintained  in  two 
characteristic  sections;  and 
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-  Some  in  plane  de-formation  is  essential  -for  hinge 
k  i  nema t  ics. 

The  model  is  given  in  Figure  5.3. 

The  rolling  radius  is  negligible  in  comparison  with  the 
other  hinge  dimensions.  Points  A,  B,  C,  and  R  are  in  a 
straight  line.   The  coordinates  of  point  B  are: 

B(x,y,z)  =  (h,   bcosp  -   [bsinp(2h-  bsinp)  ]    ,  0)       <5.4) 

From  continuity  of  the  middle   section,   the  web  height,  b,  is 
de  term  i  ned: 

/  2     2.1/2 
b  =  ZA  +  (yA  +  V  <5.5) 

With  y   =  y  ,  and  equations  <5.4),  and  <5.5): 
A     B 

z   =  bsin2p  -  hsinp  +  [bsinp (2h  -  bsinp)]    *cosp     <5.6) 
A 

An  analysis  of   the   rotation   of   the   fiber   in  the  top 

flange  will  lead  to  the  equation  for  the  hinge  length,  2h .  For 

a   fiber   which  is  initially  parallel  to  the  beam   and   passes 

through  point   A,   longitudinal   continuity   is   implied.   The 

distance   from   this   fiber   and  the  edge  KL  is  AD  =  b   -   z  . 

A 

Initially,  the   coordinates   of   the  fiber  and  the  line  GK  are 

CO,  b-z  ,  03.   After   rotation,   the   coordinates   are   [<b 
A 

z  )sin/?,  (b  -  z  )cosp,  03.  Point  A  has  the  coordinates  Ch,  y  , 
A  A  A 
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2  J.   The  distance  from  point  A   and  the  rotated  point  is  h  if 
A 

the  web  and  -flanges  are   i  nex  tensi  bl  e  .    The  equation  -for  h  is 

therefore  : 

2  2  2     2 

h  =  [h-(b-z  )sinp]   +  [y  -  (b-zA)cosp]   +  z  k  (5.7) 

Substitution  of  <5.4),  and  (5.6)  into  (5.7);  and  meeting  the 
condition  required  by  <5.5),  it  may  be  said  for  any  ratio  a/b, 
and  any  rotation  angle: 

2h  =  a   or   2h  =  b  <5.S) 

The  value  which  yields  the  lowest  energy  absorption  will  be 
used.   Thus, 


< 

2h  =  a      if  a  =  b 


2h  =  b      if  b  =  a 


(5.9) 


The   second   phase  of  deformation  is   complete   when   the 

rotation  angle  equals  the   jamming   angle,  a     .  Jamming  starts 

J 

when  the  two  buckled  halves  of  the  compression  flange   contact 

each  other.   From  Figure  5.3  the  jamming  angle  is: 


o       ,h-  0.5tx  <5.10) 

Ot   =  2arcsmp( — ; ) 

J  b 


5.1.4  Nominal  Energy  Absorbed  During  Hinge  Rotation 


The   assumptions  listed  in  the  previous  subsection   imply 
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that  all  the  energy  is  plastic  and  is  absorbed  along  the  yield 
lines  of  the  theoretical  collapse  mechanisms.  Relative 
rotation  along  the  flange  edges  JF,  and  IH  results  -from  hinge 
rotation,  see  Figure  5.1.  From  geometry,  o-f  Figure  5.3: 

TT  b 

Y  =  2  "  D  "  arcsin (1  "  ^sinp)  <  5 .  1 1 ) 

The    energy    absorbed    along    these    edges    becomes: 

W,    =  W        +  V7,  =   2m  a[—  -    p   -   arcsin  (1   -  r-sinp)  ]  <5.12> 

1  JF  IH  p      2  n 

m      =   Average    -fully    plastic    moment    of    a    unit    width    of 
P 

the    f 1 ange 
Similarly,     it    may    be    said: 


„         V7  r  _  ,,        b    .      .  ..  (5.13) 

W     =  v:        =  m  a[7T  -   2arcsin(l   -  rsxnp)  ] 
^  CR  p  h 


The       energy      absorbed    by    the    yield    1  ine    AB    can       be       found 

from    (5.11),    where    the    length    of      AB   =    z       is    Known    from    (5.6), 

A 

and  the  relative  rotation  is  ( n   -    2/J)  .  Hence: 

n  1/9  K 

w„    =  W        =  m    [bsin"p-hsinp+(bsinp  (2h-bsinp)  )         *cosp}  [7T-2arcsin  (1-r-sinp)  ] 
J  AB  p  h 

Similarly,    along      BG,       the       length       is      h,       and       the       relative 
rot  at  i  on    is    V2: 


W^    =  W„„   +  ttf        =  m  h7T  (5.15) 

4  BG  BE  p 


Figure  5.4  shows  the   relative   rotation,   r/ ,  along  yield 
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The  An^lst   ow    Rotation 


-Figure"    5,4 


line  GK.  The   construction   o-f   this   figure   is   detailed   as 

■follows:  the   plane  oo  is  normal  to  GK  and  contains  point  A. 

xy 

Point  A'   is   the   projection   o-f   A  onto  the  xy  plane.   A   is 

o 

obtained  by  rotating  to  abou  t  its  intersection  co  with  the  xy 

xy 

plane.       Coordinates    o-f    A"    are: 

1/2 
htanp  +  bcosp   -    [bsinp (2h-bsino) ] 

y   „    =  ,  2  <5.16) 

*A"  1   +   tan    p 

V   =  yA»tanp  (5.17) 

The    energy    absorbed    along    both    the    GK,    and    EL    1  ines    are    then: 

_£ i    (5.18) 


W5  =  V  +  WEL  "  2»  barctanl— —  , 

[(h-xA„)  +  (yA„-yB)-] 

The  rolling  radius  along   the   yield   lines   GA,   and   AE 
remains  relatively  constant,  so: 


w*  =  w™  +  w.^  =  2m  (-)z  (5.19) 

6     GA     AE      p  r   A 


During   hinge   rotation,   the  rolling   radius   varies   as 
descr  i  bed  by : 

r  =  (0.07  -^)h  <5'20) 

Equation  (5.20)  is  based  on  empirical  data.  Kecman 
approximated  the  variation  o-f  the  rolling  radius  along  KA  with 
a  hyperbolic  approximation: 

r   =  «  (5.21) 

KA    £K 
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Where  1   is  the  distance  from  K  along  KA.  The  linear  variation 
k 

o-f  the  rolled  length  between  K  and  A  is  given  by: 


JS  z  (5.22) 

~r    KA  A 

The  energy  along  KA  is: 

2m  z^KA 
W   =    P  A —  <5.23) 

KA       3r 

or 

W7  =  wka  +  wta  =  Tm  — fh2  +  y  2  +  z  2]1/2  <5.24) 

'     KA     LA    3   p  r       -^B      A 

The  relative  rotation  along  yield  lines  PN ,  and   QM  \sp, 
whereas  along  KL,  and  MN  the  angle  is  given  by: 

t,  =  arctan  ( — ) 

yA 

Finally,       the       energy         absorbed         by       these       yield       lines       is 
approximated    by: 

W8   =   %+WQM+WKL+WMN   =    2mp[ap   +   2harctan  (%  ]  <5.25> 

'A 

At  this  point,  the  total  energy  absorbed  by  the  hinge  may 
be  expressed  as: 

I        ,   ,  <5.26) 

w(a)    =  Z  W.  (a) 

i=l 
This    result     is    valid    for    a      rotational       angle       less       than       the 

jamming    angle.       Kecman    determines    the    moment-angle 
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relationship    numerically    -from    the    -following    equation: 

M(a)    =  W<Q  +  Aa)   -  w(a)  <5.27) 

Where  -A<^  represents  a  small  but  -finite  increment   o-f  the  hinge 
angl e ,  a. 

This  -formulation  assumes  that  the  material  is  ductile 
enough  not  to  -fracture  during  hinge  rotation  or  subsequent 
bending  collapse.  This  is  considered  a  sa-fe  assumption  -for 
common  structural  steels.  The  average  plastic  moment  is  given 
by  the  -following  equation: 

a  t2 
mD  =  ou  <5.28> 

p      4 

The  yield  stress,  and   the   maximum   nominal   flow   stress   as 
-functions  o-f  temperature  were  determined  in  Chapter  5. 

5.1.5  The  Final  Form  o-f  the  Hinge  Moment 

The  analysis  in  the  preceding  section  is   not   valid   -for 

the   -first   phase  oi         de-formation.     This    phase   will   be 

approximated  by  assuming  that  -for  a  moment   greater   than   the 

•fully  plastic   moment,   M  ,   is  constant  and  equal  to  M  .  This 

o  o 

would  be   a  conservative  assumption.   For  hinge  angles  greater 

than  the  jamming  angle,  Kecman  o-f  tens  the  following  -formula: 


M(a)  =  M(a  )  +  1.4  [M  -  MCoij)  ]  (a-  Oj) 


(5.2?) 
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This  t'ormu  1  at  i  on  has  been  programmed  in  a  M  i  crosof  t-BASI  C 
program.  The  source  code  is  given  in  Appendix  B.  Figure  5.5 
shows  the  results  of  the  program. 


5^2_tiaxduk.Ls— and_U±£czb±ck±.l3_±l£.±h£id 

The   primary  difference  between  Kecman's  method  and   this 
method   is   that   Hayduk   and  Wierzbicki  considered   not   only 
i  nex  tensi  onal    de-formation    modes,   but    also    ex  tensi  onal 
de-formation   patterns.   The  e-f-fect  o-f  these  extensional   modes 
will  be  -fully  developed  in  this  section. 

Hayduk's  and  Wi erzbi ck i ' s  -formulation  results  in  a  closed 
form  analytical  solution  which  predicts  the  mean  crushing 
force.  Reference  11  specifically  addresses  the  problem  of 
analyzing  the  crushing  process  of  thin  walled,  plate  formed, 
open  structures  with  particular  emphasis  on  L  shapes  and 
cruci forms.  The  I  beam  can  be  approximated  by  two  T  shapes  or 
four  L  shapes.  Lower  and  upper  bound  solutions  for  the  mean 
crushing  strength  of  cruciforms  was  obtained  by  considering 
various  modes  of  deformation.  According  to  Hayduk  and 
Uierzbicki,  extensional  deformation  accounts  for  at  least  one 
third  of  the  dissipated  energy. 
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The  end  products  of  this  section  will  include  the 
application  o-f  Hayduk's  and  Uli  erzbi  ck  i  '  s  method  to  an  I  beam, 
and  reviews  of  both  references  11  and  12. 


5.2.1  Collapse  Modes 

Hayduk  and  Wierzbicki  considered  out  o-f  plane 
deformations  in  which  the  deflections  exceed  10-100  times  the 
section  thickness  and  are  of  the  order  of  magnitude  of  the  web 
height  and  the  flange  width.  The  process  of  collapse  starts 
as  the  compression  -flange   loses  stability,  either  elastically 

or  plastically  at  a  critical  load,  P  .  Following  the   buckling 

c 

deformation,   the   flange   can   maintain  added  loads   until   a 

maximum  load,  P    ,  is  reached.   Eventually,  as  the  resistance 

max 

decreases,  the  structure  collapses  completely.   The  structures 

being  examined  here  often   regain   strength   after  the  initial 

loss    of   resistance.    Subsequent   load   peaks   reveal    the 

development  of  additional  folds. 
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5.2.1.1  Isometric  Trans-formation  Modes  <  Inex  tensi  onal 
Def ormat  i  ons) 

These  modes  are  of  the  type  used  by  Kecman  to  analyze  the 
collapse  mechanism.  In  this  type  of  deformation,  the  energy 
is  concentrated  over  hinge  lines,  while  the  rest  of  the 
structure  remains  rigid.  An  assumption  of  i nex tensi onal 
deformation  is  that  material  continuity  is  maintained.  During 
the  transformation  of  two  surfaces,  the  lengths  of 
corresponding  arc  elements  is  preserved.  Hayduk  and 
Ulierzbicki  show  the  progression  of  various  deformation  paths. 
These  are  illustrated  in  Figure  5.6. 

In  Figure  5.6  <a),  and  <b),  the  final  state  is  formed 
i nex tens i bl y .  This  can  only  be  done  if  the  boundaries  have 
sufficient  freedom  to  deform.  The  element  in  Figure  5.6  <c) 
is  initially  folded  along  the  stationary  hinge  line  AB ; 
finally  the  element  is  rotated  extensibly,  <ref.  11). 


5.2.1.2  Hinge  Collapse  Modes 

There  are  two  types  of  hinge  collapse  modes.  The  first 
is  a  propagating  hinge  1  i ne  which  occurs  when  the  hinge  1  i ne 
is  free  to   move   with   respect   to   material   points,  and  the 
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material  o-f  one  -flange  is  ef  -fee  t  i  ve  1  y  transmitted  to  another 
■flange.  Extensional  de-formation  occurs  in  the  vicinity  o-f  the 
hinge  line.  This  mode  is  a  characteristic  o-f  thin  walled 
structures  consisting  o-f  two  intersecting  plates  such  as  the  I 
beam  being  considered  in  this  study. 

The  second  hinge  collapse  mode  is  the  formation  o-f  a 
stationary  hinge  line.  This  mode  is  characterized  by  a  hinge 
line  fixed  in  the  material.  Side  flanges  which  are  initially 
rectangular  are  transformed  into  trapazoidal  elements. 
Therefore,  extension  occurs  over  an  entire  flange.  This  mode 
will  not  be  considered  in  this  analysis,  <ref.  11). 


5.2.1.3  Assembled  Collapse  Modes 

Hayduk  and  Wierzbicki  demonstrated  the  method  of  assembling 
the  folding  patterns  shown  in  Figure  5.6.  The  patterns  are 
folded  into  more  complicated  symmetric,  asymmetric  and  mixed 
collapse  modes  which  are  shown  in  Figure  5.7. 
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5.2.2  Anal ysi  s 

For  the  purposes  o-f  this  stud/,  the  analysis  o-f  a  T 
section  will  be  used;  this  analysis  is  considered  to  be 
equivalent  to  an  I  beam  analysis.  In  re-ference  12,  the  short 
T  cross  section  under  consideration  was  subjected  to 
compression  with  clamped  boundary  conditions.  A  model  is 
presented  in  Figure  5.8. 

For  the  purpose  o-f   estimating   the   mean  crushing  -Force, 

P  ,  the  T  section  is  expanded  into  a  plate   and   an   angle   as 
m 

shown  in  Figure  5.9. 

The  plate  will  collapse  through  the  bending  zone  limited 
by  the  two  horizontal  moving  lines.  The  angle  will  -follow  the 
-folding  mechanisms  described  by  Figure  5.10(b).  The  collapse 
modes  in  Figure  5.10<b)  is  consistent  geometrically  and 
k i nemat i cal 1 y .   The  regions  are  described  as  follows: 

(I)  Four  plane  trapazoidal  elements  as  rigid  bodies. 

(II)  Two  sections  o-f  cylindrical  surfaces  at  which 
continuous  bending  takes  place  without  any  extension. 

(III)  Two  sections  o-f  conical  surfaces  in  which  material 
is  bent  and  rebent  again  as  the  material  moves  -from  one 
■flange  to  another. 
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<IV)  A  section  of  toroidal  sur-face  which  produces 

extension  in  a  c  i  rcumf  ren  t  i  al  direction  and  continuously 
changes  principal  curvature  in  the  other  direction,  <ref  . 
11) 

The  -following  assumptions  will  be  used  in  this  analysis: 

-  Neglect  the  short  side  o-f  the  angle  element. 

-  The  length  o-f  the  two  horizontal  hinge  lines  is 
approximately  equal  to  the  length  o-f  the  large  side  of 
the  angl e  e 1 emen t . 

-  Two  incl  ined  hinge  1  ines  are  of  a  length  equal  to  a 
hal-f  wave  length,  H. 

-  The  length  of  the  toroidal  surface  in  the  vertical 
direction  is  neglected. 

The  energy  dissipated  of  an  angle   element   is   given   by 
-ummation  of  three  terms,  <ref  11): 

1.  The  energy  of  a  continuous  deformation  field  in  the 
section  of  a  toroidal  shell: 

El  -   9-28MoHr/tw  < 5  .  30 ) 

2.  The  energy  dissipated  by  the  discontinuous  velocity 
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field    in    the    horizontal    hinge    line    of    length,     1: 

E2=§V  C5.31) 


3.  The  work  done  by  an  incl  ined  hinge  1  ine  through  the 

de-formation  process: 

2 
E,  =  2.22M  H  /r  •  =  o-^ 


where  : 


r  =  The  small  radius  of  the  toroidal  surface. 

2 
M   =  <^Ct  ]  /4,  the  fully  plastic  yield  moment  per  unit 
o     o   w 

width 


H  =  Mode  half-wave  length 


a      =  0.7a 
o        ou 

a         =    Ultimate    strength 
ou 

t      =   Web    thickness 
w 


For    the    models    given    here: 

2HP      =   E,    +   4(E_)       .     +   4(E    )  ,.  +    2E 

m  1  2  web  2   flange  3  (5.J3) 

After    substitution    of    the    dissipated    energies    into    (5.33): 

P  /M     =   4.64f     +  2L[C   +  d(-£)    ]    +   2.22-^  <5'34> 

\        H  tw  r 


-    77    - 


Optimizing    equation       <5.34)       by      using       the       -following    partial 
der  i  vat  i  ves : 


3P 
a   m 


=    0 


(5.35a) 


d   m 

3h 


=   0 


Equating    (5.35a)    and    (5.35b)    yields: 


r2.22  2,  *f    2)]1/3 

r   =    [ irt       (c   +  d  — )     ;  J 

4 .  64  w 


2 

7T    t 


H   =       [T 


V2.  2.1/3 


w f   z   2 

"(2.22)  (4.64)'    (c   +  d  (— >     )    3 


(5.35b) 


(5.36) 


(5.3?) 


w 

Ultimately,  the  three  basic  mechanisms  o-f  energy 
dissipation  appear  equally  weighted  in  the  -final  equation  -for 
mean    crushing    -force: 


t     2 

P   /M      =    3[(4.64)  (2.22)7T(1    +   k  (— )     )  (£-)  ]    / 
mo  t  t 

w  w 


(5.38) 


Equation    (5.38)    may    be    simplified    as    -follows 


V  Mo   "    3 [A1A2A3  ^ ] 1/3 

w 


(5.3?) 


where  : 


k  = 


(5.40) 


A2    =    4.64 


(5.41) 


V2 


A2    =   Tfl    +   k(— )     ] 


w 


A3    =    2.22 


(5.42) 


(5.43) 


Note    also    that    equations    (5.34),     (5.36),       and       (5.37)       may      be 
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simplified    into    the    t'oll  owing    forms: 


P  /M      =   A, +   A  —  +   A  —  f«=i    44) 

m      o  It  2H  3r  l3'^' 

w 


2    3      1  w 

,2,.    ...   2  1/3  (5.46) 


"-    HA^A^Xc^)] 


The  appl  i cat  ion  of  these  equations  to  the  case  of  an  I 
beam  requires  that  the  geometric  parameters  be  modified.  An  I 
beam  is  two  T  sections,  hence  the  dimensions  are  doubled;  but, 
b   =    2c,    which    yields    the    following    relations. 

VMo   =    3  [A1A2A3  (b/2V  ]  V3  <5.47) 


k=2a/b  <5'4S> 


It  is  now  possible  to  apply  this  relation  determined  for 
the  crushing  case  to  the  bending  case.  The  geometry  of  this 
crushing  mode  is  given  in  Figures  5.11  and  5.12. 

Recall  the  propagating  hinge  line  mode  as  described  in 
paragraph  5.2.1.2.  This  mode  describes  the  global  geometry  of 
the  crushing  process  and  is  determined  purely  by  the  width  of 
the  two  side  panels,  d,  which  is  prescribed,  the  wavelength, 
H,  which  is  to  be  determined,  and  the  crushing  rotation  angle, 
or,  which  is  a  parameter  of  the  process.  The  relation  between 
the  axial  displacement,  <*>  ,  and  the  angle,  <*,  is: 

6  =  2H(1  -  cos6)  <5.4?) 
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p 


p 


Crushing  £asi= 


V\w$e     5".U 


ExTCNSlONAU    CaLLAPSe   V\oi>£  - 

^fiHmoN  op=  Geometric  f^RA^eTERS 


fn^uae     5.IZ. 


Re -fere  nee  12  determined   an   expression   -for  the  instantaneous 

crushing  force,  P,  normalized  with  respect  to  P  : 

m 

lr„        o  2  0.9cosB  ,  <5.50) 

P/Pm  =  ^cosB  +  ^r^  +  tan3(1   +   sin^B)i/^ 


Figure    5.13    gives    the    geometry    of    the    bending    case.       From 
Figure    5.13: 


.      a       6/2 

tan2    =— 


(5.51) 


Thus,  it  may  be  said  that: 


P/P  =  f(6)  =  f(2btar£)  (5.52) 

m  2. 


For  a  section  width,  b: 


M  =  -  b  (5.53) 

2 

P 
M  =  —  b  (5.54) 


m 


n 


Where  M   is   the   mean   moment   due  to  bending  rotation,  a,    as 
m 

de-fined  in  Figure  5.14.  After  equating  (5.53),  and  (5.54): 

§-  =  £-  (5.55) 

M     P 
m    m 


Figure       5.15      gives    P/P      as    a    -function    o-f    <V2H      based      on 

m 

equations    (5.49),    and    (5.50),     (re-f.       12).       This      graph       is      an 

equivalent         description         o+"         M/'M  as         a         -function         of 

m 

(b/H) C  tan(a/2)]  . 
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Figure  5.  !5" 


For  this   purposes  o-f  this  study,  an  approximate  relation 

t'or  M/M   and  a    will  be  determined.   Assume   equation  (5.56)  is 
m 

a  valid  -function  -for  M/M  : 

m 


M/M   =  D/(6/2H)n 
m 


(5.56) 


where  : 


D,  n  are  constants 


Assume  the  new  curve  goes  through  points  P  ,  and  P  ,  in  Figure 

1        2 

5.15. 


P.  (5/2K,  M/M  )  =  (0.30,  1.0) 

x  m 


<5.57) 


P_(6/2H,  M/M  )  =  (1..0,  0.24) 
z  m 


•:5.53) 


Sol  m  i  ng  -for  D,  n  : 


D  =  0.24 


(5.5?) 


n  =  1.2 


(5.60) 


Assuming    n    =    1.0,     the    -following    equation    results: 


M/M       = 


0.24 


(5.61 ) 


m  "     6/2H 

Figure  5.16  shows  how  the  approximate  solution  compares  to  the 
exact  solution.  Substituting  equation  (5.51)  into  equation 
(5.61),  the  following  results: 
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z  - 


m 


I    -■ 


Z^r 


■  is 


l.o 
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o  24  ( 5  • 62) 

M/M      = 


m        (b/H)tan(a/2) 
Considering    equation    (5.54): 

Pmb   0.24 (5.63) 

2     l(b/H)tan(ct/2)] 

Now  substitute  equation  (5.47)  into  (5.63) 

M=    [3M    [A.A.,A_(b/2t    )]1/3]  [b/2]  [°'24    H ]  (5.64) 

°      1    2    3  w  '    J  lbtan(a/2)  J 

But    M       is    the    -fully    plastic    yield   moment    -for    the    flange 
o 

2 

ao  (5.65) 

M  = 


o  4 


where    t    =    t      =    t     . 
-f  w 


The  desired  -final  result  o-f  this  analysis  is   a   relation 

between  M  and  M  ' .  Rewriting  equation  (5.65): 

o 


M  '  =  a  j    zdA 

o      o  A 


(5.66) 


Recall     equations    (2.2)    and    (2.4): 

a    (T)    =    -20. 5T    f    40,000.0  (2.2) 


o 

J  zdA  =  abt  +  ?-L 

A  4 

Finally,  substitute  equation  (2.4),  (5.66),   (5.65),   and 
(5.46)  into  (5.64) : 
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I 

o 


J 


<r 


\x/h£R2.gickis  Analysis 


2  ?  2         1  /? 

30    t  0.24[(A!/A   A    )b"2t]i/J 

M/M°'%0o(abt,b^,'AlA^(b/2t"      '|H      *-»"> ' 

The    simplified    final     result    of    this    relation     is: 

,         0.36TTt  (b  +   2a)  (5.68) 

V   o      ~~   tan (a/2) [b(b  +  4a)] 

Figure  5.17  shows  graphically  the  resulting  moment-angle 
rotation  relationship  as  determined  by  Hayduk  and  Ulierzbicki. 
Appendix  C  contains  the  source  code  -for  a  Microsoft  BASIC 
program  which  determines  the  moments  resulting  from  equation 
<5.72>  . 


■5^3-L.Qmp  a  c_l.s,q.q— of — &£=,u±±=,-£>£-.±h£.-.  Jljjd— tlelbcud-s 

The  primary  difference  between  the  two  methods  of 
analyses  is  that  Kecman  considered  only  i nex tens i onal 
deformation  modes,  whereas  Hayduk  and  Ulierzbicki  considered 
extensional  deformations  as  well  as  i nex tens i onal 
deformation.  Since  the  second  method  considered  more  energy 
dissipation  than  the  first,  the  allowable  moments  due  to 
deformation  are  greater.  However,  for  angles  greater  than  one 
tenth  of  the  jamming  angle,  Kecman 's  method  yields  greater 
moments.   The   errors   induced   by   his   numerical  integration 
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technique  are  compounded  as  the  angles  increase.  Figure  5. IS 
shows  graphically  how  the  results  of  each  analysis  differs. 
The  dimensions  of  the  section  analyzed  are  given  below: 

a  =  1 2  i  nches 

b  =  1 2  i  nches 


t   =  t   =  . 6    i nches 
f     w 

The  maximum  moment  given  in   Figure   5.13  is  the  fully  plastic 

yield  moment  of  the  section,  hereafter   called  M  ,   which   is 

o 

only  a  function  of  geometry  and  temperature  of  exposure  of  the 

beam.   Hence,  M   is  not  a  function  of  the  analysis  method. 

o 


Clearly,  as  the  angle  of  rotation  increases  the  moment 
decreases;  and  as  the  temperature  increases,  the  moment 
decreases.  Appendices  B  and  C  also  contain  results  which  show 
the  effects  of  changing  the  geometry  of  the  section.  For  both 
analyses,  the  moments  are  decreased  as  the  section  thickness 
decreases;  the  web  height  decreases;  or  the  flange  width 
decreases . 

These  methods  of  analyses  are  considered  accurate 
predictions  of  the  failure  mechanisms  of  a  beam  exposed  to 
fire.  Recommendations  for  using  the  information  provided  by 
these  analyses  will  be  given  in  the  following  chapter. 
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Chapter    6 
rnhiri  ms tons  akjh  prrnMMFNnAT t dns 


A  summary  of  each  chapter  is  presented  here  with 
concluding  remarks.  Finally,  recommendations  are  offered  for 
using  this  analysis  -for  -future  studies. 


As  stated  in  the  introduction,  the  ultimate  objective  o-f 
this  study  was  to  show  how  much  the  strength  of  a  shipboard 
structural  member  is  diminished  when  exposed  to  fire.  To  this 
end,  Chapters  2,  3,  and  4  provided  input  into  the  analyses  of 
the  failure  mechanisms  studied  in  Chapter  5.  The  results  of 
the  preceding  chapters  will  be  summarized  here. 

Chapter  2  developed  the  model  of  the  structure.  Simple 
observations  revealed  that  the  ship  is  a  redundant  structure 
composed  of  beams,  plates,  and  columns.  The  motivation 
towards  choosing  an  I  beam  for   the   analysis   was   developed. 
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The  material  selected  -For  the  model  was  ASTM-A36  steel  as  its 
physical  and  mechanical  properties  are  similar  to  the 
properties  of  ABS  Grade  A  steel,  a  common  shipbuilding  steel. 
The  equations  -for  the  temperature  dependent  mechanical 
properties  were  the  output  of  this  chapter  to  be  used  in  the 
anal  ys i  s . 

Chapter    3    provided   background    for    characterizing 

shipboard  fires.   Three  parameters  which   are   assumed   to   be 

accurate  predictions  of  the   contributing   factors   of  a  fully 

developed  fire   were  defined.   These  factors  are:  the  duration 

of  the  fully  developed  fire;  the  overall  penetration  flux;  and 

the  average  fire  temperature.   Since  fire  is  indeed  a   complex 

process  with   several   unquan t i f i abl e    variables,    a    time 

temperature  curve   was   used   as   a   definition   of   the   fire 

process.   For  the  purposes  of   this   study,  a  step  function  to 

o 
1700  F  was  used  as  the  time  temperature  relation. 

Chapter  4  was  the  heat  transfer  analysis  which  determined 
the  relationship  between  the  elevated  temperature  exposure  and 
the  distance  along  the  length  of  the  beam.  The  results  of 
this  analysis  led  to  the  conclusion  that  a  beam  exposed  to 
fire  will  experience  the  same  failure  mechanisms  as  a  beam 
under  a  point  load. 

Finally,  the  primary  thrust  of  this  study  was  to   analyze 
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the  failure  mechanisms  due  to  -fire  exposure.  Two  methods  were 
•fully  developed  in  the  preceding  chapter  and  the  results  were 
summarized  graphically  in  Figure  5.13. 


*-?    Pnrnmmflnfjftti,Qng 

The  ultimate  objective  o-f  any  structural  design  process 
is  to  develop  a  structure  which  is  strong  enough  to  withstand 
all  expected  service  loads  and  a  high  percentage  of  unexpected 
loads.  The  information  provided  by  the  analysis  in  Chapter  5 
is  considered  useful  input  for  the  design  of  shipboard 
structures.  It  is  possible  through  the  use  of  Figure  5.18  to 
determine  how  the  stiffness  of  a  structural  member  is 
decreased  when  exposed  to  elevated  temperatures.  The 
diminished  stiffnesses  may  be  used  as  input  to  a  finite 
element  program  which  simulates  a  ship  structure.  The  effect 
of  a  shipboard  fire  upon  the  entire  structure  may  then  be 
de  term  i  ned . 

The  examination  of  the  effect  that  one  weakened  member 
has  on  the  other  members  of  a  redundant  structure  is 
considered  fertile  ground  for  future  studies.  Reference  16 
may  provide  insight  into  the  behavior  of  this  phenomenon. 
Zubelewicz  and  Mroz  numerically  simulated  rock  burst  processes 
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treated  as  problems  of  dynamic  instability.  According  to 
reference  16,  local  plastic  failure  occurs  at  one  location 
(such  as  in  a  member  of  a  ship's  structure)  and  the  elastic 
energy  flows  to  the  hinge  and  is  the  driving  force  into  the 
damage.  Energy  then  flows  into  the  neighboring  hinges,  which 
is  to  say  that  one  unstable  member  may  produce  a  chain 
reaction  and  other  members  may  subsequently  fail. 

As  a  final  word,  it  is  highly  recommended  to  utilize  the 
information  provided  in  this  study  as  an  aid  for  examining  the 
effect  of  fires  on  the  ship's  structure.  Ultimately  a  ship 
may  be  designed  to  withstand  shipboard  fires. 
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Appendix  A 
EE^I-IEANSEEE-AN^LXSlS-EEOGEAfcl 


10  REM  PROGRAM  TO  CALCULATE  TEMPERATURE  GRADIENT  IN  A  BEAM 

20  DIM  TEMP<20) ,Y(20) , Z < 20  ) ,TGP<20) 

25  INPUT  "ENTER  THE  BEAM  LENGTH  IN  FEET" ;XL 

30  INPUT  "IS  BEAM  SOLID  (S)  OR  I  ( I ) " ;N* 

40  IF  N*="I"  THEN  GOTO  90 

45  LPRINT  '  SOLID  BEAM" 

50  INPUT  "ENTER  WIDTH  OF  BEAM" ;W 

60  INPUT  "ENTER  DEPTH  OF  BEAM";D 

65  LPRINT  "  WIDTH,  DEPTH=" ;  W;D 

70  B=SGR< ( 12/5.4>*<<l/D)+< 1/W))) 

80  GOTO  155 

84  LPRINT  "  I  BEAM" 

90  INPUT  "ENTER  WIDTH  OF  FLANGE" ;W 

100  INPUT  "ENTER  DEPTH  OF  WEB";D 

105  LPRINT  "  WIDTH  OF  FLANGE,  DEPTH  OF  WEB=" ;  W;D 

110  INPUT  "ENTER  THICKNESS  OF  FLANGE" ;T1 

120  INPUT  "ENTER  THICKNESS  OF  WEB" ;T2 

125  LPRINT  "  THICKNESS  OF  FLANGE,  THICKNESS  OF  WEB=" ;T1 ;T2 

130  P=2*<D+2*W+2*T1-T2) 

140  S=2*(T1*W)+(D*T2) 

150  B=SQR< <12/5.4>*<P/S>) 

155  LPRINT  "  LENGTH  OF  BEAM=";XL;"  FEET" 

156  LPRINT  "  B=";B  :  LPRINT  "  " 
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160    BL=B*XL 

170    BOT=EXP<BL)-EXP<-BL) 

180    X=XL/20 

190  LPRINT  "  PT  NO  X  T-Ta/Tb-Ta" 

200  FOR  KX=1  TO  20 

210  Y<KX)=X*KX 

220  Z<KX)=B*<XL-Y<K7.>) 

230  TOP<KX)=EXP<Z(KX))-EXP<-Z(KX)) 

240  TEMP<KX>=T0P<KX)/60T 

250  LPRINT  "  M ;KX,Y<K%) ,TEMP<KX> 

260  NEXT  KX 

270  END 
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Appendi  x    B 


KE£mkHS-.&U&LX3±S-£EDBE&tl 


10    DIM   W<20) 

20  INPUT  "TEMPERATURE  OF  EXPOSURE,  DEGREES  F" ;T 

30  INPUT  "FLANGE  WIDTH,  INCHES" ;A 

40  INPUT  "WEB  HEIGHT,  INCHES" ;B 

41  INPUT  "THICKNESS,  INCHES" ;TH 

42  LPRINT  "  " 

50  LPRINT  "  " 

51  LPRINT  "  TEMPERATURE=" ;T 

52  LPRINT  "  FLANGE  WIDTH=";A 

53  LPRINT  "  WEB  HEIGHT=" ;B 

54  LPRINT  "  THICKNESS*" ;TH 
60  S0=-20.5*T+40000 .0 

65  IF  T>=1?51  THEN  S0=4.5 

67  M0=  S0*B*TH*<A+B/4) 

68  LPRINT  "  FULLY  PLASTIC  MOMENT,  INCH-POUNDS" ; MO 
70  J=0 

90  U=3. 141 59265 
100  IF  A  <=  B  THEN  H=A/2 
110  I F  A  >  B  THEN  H=B/2 
120  Ql=<H-.5*TH>/6 

130  D=2*<ATN(Q1/SQR<-Q1*Q1  +  1  ) >> 

131  LPRINT  "  JAMMING  ANGLE,  RADIANS" ;D 

132  LPRINT  "  ANGLE;  MOMENT/MO;  MOMENT" 
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140  X=SU*<THA2)/4 

150  FOR  K=D/10  TO  11*D/10  STEP  D/10 

151  PI  =  K/2-K's3/48+K"5/3340 

152  P2  =  1 -KA 2/8+ KA 4/334 

153  P3  =  B*P1 

154  P4  =  <2*H  -  P3) 

155  P5  =  B*P2 

156  P6  -  H*P1 

157  P?   =  P1/P2 

165  YB  =  P5-<SQR<P3*P4>) 

170  R=(  .0  7-K/70)*H 

180  ZA=(P3*P1 )-P6+(SQR(P3*P4))*P2 

190  Q1=1-<B/H)*P1 

200  Q2=ATN<Q1/SQR<-Q1*Q1+1 )) 

210  Wl=2*A*<U/2-K/2-Q2) 

220  W2=A*<U-2*Q2> 

230  W3=<P3*P1-P6+<SQR<P3*P4> ) *P2> *< U-2*Q2) 

240  W4=H*U 

250  YC=<H*P7+P5-<SQR<P3*P4) ) )/< 1+P7A2) 

260  XC=YC*P7 

270  Q1=ZA/<SQR< ( ( H-XC) A2> + < ( YC-YB) A2) ) ) 

230  Q3=ATN<Q1) 

290  W5=2*B*Q3 

300    UJ7=4/3*<ZA/R)*<SQR<  <  HA2)  +  ( YB'S2)  +  <  ZAA2)  )  ) 

310    UI6=2*(H/R)*ZA 

320    Q1=ZA/YB 

330    Q3=ATN<Q1) 
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340  WS=2*<A*<K/2>+2*H*Q3) 

350  J=J+1 

360  UK  J)=W1  +W2+W3+UI4+W5+W6+W7+W8 

370  NEXT  K 

380  FOR  I7.=  l  TO  J-l 

390  Y=<TH/<B*<4«A+B))>*<W<I%+1  >-W<  I  JO  )/<D/10) 

391  M1=Y*M0 

400  LPRINT  "  " ;IX*D/10;Y;M1 

401  NEXT  17. 

402  LPRINT  "  " 

403  LPRINT  "  " 

404  LPRINT  "  " 
410  LPRINT  "  " 
420  END 
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Appendix    C 
U&XnUK-&UD-hl±ERZB±£K±±S-.&hl&LXS±S-.ERQQE&£l 


10  INPUT  "TEMPERATURE  OF  EXPOSURE,  DEGREES  F" ;T 
30  INPUT  "FLANGE  WIDTH,  INCHES" ;A 
40  INPUT  "WEB  HEIGHT,  INCHES" ;B 
50  INPUT  "THICKNESS,  INCHES" ;TH 

61  LPRINT  "  " 

62  LPRINT  "  " 

70  LPRINT  "  TEMPERATURE  " ;T 
SO  LPRINT  "  FLANGE  WIDTH" ;A 
90  LPRINT  "  WEB  HEIGHT" ;B 
100  LPRINT  "  THICKNESS" ;TH 
120  Tl=<A*B*TH+B*B*TH/4> 
130  T2=-20.5*T+40000 
135  IF  T>=1?51  THEN  T2=4 . 5 

137  M0=T2*T1 

138  LPRINT  "  FULLY  PLASTIC  MOMENT,  INCH-POUNDS"  , -MO 
140  T3=0. 36*3. 141 593*TH*< B+2*A> 

150  T4=1/<B*<B+4*A>) 

170  LPRINT"  ANGLE;  MOMENT/MO;  MOMENT" 

171  IF  A<=6  THEN  J=A/2 

172  IF  A>B  THEN  J=B/2 

175  Q1=<J-.5*TH)/B 

176  D=2*(ATN<Q1/SGR<-Q1*Q1+1 ) )) 
130  FOR  K=D/10  TO  D  STEP  D/10 

-  104  - 


131  T5=TAN<K/2) 
190  M=T3*T4/T5 
195  M1=M*M0 

200  LPRINT  "  " ;K;M;M1 

201  NEXT  K 

202  LPRINT  "  " 

203  LPRINT  "  " 

204  LPRINT  "  " 
210  LPRINT  "  " 
220  END 
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